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We discuss CP violation in top-antitop production at the LHC, induced by gluon 
fusion and hnal-state Higgs exchange. Results by Bernreuther and Brandenburg are 
confirmed and further reduced. The lepton energy asymmetry is studied in detail in 
explicit Two-Higgs-Doublet Models with near-maximal mixing in the neutral Higgs 
sector. Unless there is only one light Higgs particle, and unless (in "Model H") 
tan/3 ^ 1, the CP-violating effects are very small, possibly too small to be seen at 
the LHC. 

1 Introduction 

One of the most promising ways in which one can search for new physics at the LHC, is in 
CP violation in connection with ti production. It is generally believed that the top quark, 
since it is very heavy, might be more susceptible to new physics [1,2]. In the particular 
case of Higgs-mediated interactions, this is naturally the case, since the Higgs coupling to 
the top quark is proportional to its mass. This process 
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has been explored in considerable detail by Bernreuther and Brandenburg [3,4] who iden- 
tified the different kinematical structures appearing in the CP- violating part of the inter- 
action, and evaluated them in a generic Two-Higgs-Doublet Model [5]. 

Here, wc review (and confirm) these calculations, and apply the results to a particular 
version of the Two-Higgs-Doublet Model, in which the CP violation is minimal in structure 
[6]. This allows us to relate and constrain the couphngs and masses of the model. Such 
relations among parameters are crucial in order to estimate the magnitude of possible 
signals. 

If the Higgs states are not eigenstates under parity, then their couplings to the fermions 
will violate CP. In particular, if the Hti coupling (for a given Higgs particle) is of the 
generic form 

Hti: [a + i75a], (1.2) 

then the CP-violating part of the cross section for the process (1.1), which depends on the 
t and i spins, will be proportional to the dimensionless model-dependent quantity 

TCP = -oo (1-3) 

where a, a are the reduced scalar and pseudoscalar Yukawa couplings, respectively. 

Here, we restrict ourselves to the subprocess gg — > tt since it is the leading ti production 
mechanism at the CERN LHC. The qq initial states, which are important at the Tevatron 
[3], contribute at the LHC less than 10% to the total cross section [7] and give only a 
numerically negligible contamination of the CP- violating signal [1,3]. 

The MSSM provides an alternative, very interesting framework for CP violation via 
the Higgs sector [8]. Additional effects in the MSSM include CP- violating gluino exchange 
[9, 10] as well as those due to phases in the bilinear and trilinear couplings, given by the 
so-called ji and At parameters. For an apphcation to the process (1.1), see [10]. We here 
restrict ourselves to the 2HDM, this is already a rich framework. 

The paper is organized as follows. After a review of notation and relevant formulas in 
sect. 2, we give model-independent results in sect. 3, focusing on basic observables discussed 
in [3]. In sect. 4, wc review the physical content of the 2HDM, and in sect. 5, we study 
the magnitude of the CP violation for two distinct neutral Higgs mass spectra: two light 
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and one heavy vs. one light and two heavy. Sect. 6 is devoted to concluding remarks, and 
an appendix summarizes the basic one- loop results for the CP- violating amplitudes. 

2 Preliminaries 

A schematic representation of a generic one-loop diagram of the process gg —>■ tt is given 
in fig. 1.^ 




t 



Figure 1: Kinematics of the underlying g{ki) + g{k2) — » + ^(^2) reaction. 
2.1 Notation 

The momentum and the spin four-vectors of the top (and antitop) quark are denoted by 
Pi and Si (and p2 and S2), respectively, with nit the mass, whereas the momenta and the 
Lorentz indices of the initial gluons are represented by ki and /ij {i = 1,2). We introduce 
the linearly-independent set of momenta {Pg, Q, Pf): 

Q = ki + k2=pi+P2, Pg = ki-k2, Pt=Pi-P2, (2.1) 

with Q ■ Pt = Q ■ Pg = 0. The non- vanishing scalar products involving the momenta {Pg, 
Q, Pt) are given by = —Pg = s, Pt = —[3'^s and {Pg ■ Pt) = —jdzs where s is the 
gluon-gluon center of mass energy squared, (3 = a/1 — 4m^/s is the top quark velocity 
and z = cos 9 = (Pg ■ Pt), with 9 the scattering angle in the gluon-gluon center of mass 
frame. Working with the linearly-independent set of momenta {Pg, Q, Pt) simplifies the 
kinematics, the four-point first- and second-rank tensor loop integrals in particular. 
^To produce the figures, the AXODRAW package [11] was used. 
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Figure 2: Lowest-order QCD Feynman diagrams of the underlying gg — > tt reaction, the 
crossed diagram of (a) is not shown. 



Denoting by M.q the QCD amphtude corresponding to the diagrams of fig. 2, and by 
M-i the one corresponding to the one-loop Higgs-exchange diagrams of fig. 3, we obtain 
the cross section 

croc |Mo + Mil'- (2.2) 

Spin-dependent CP- violating parts arise from terms linear in 7cp, (1.3). Those originate 
from interference between Jv[q and A^i, and from |A^ip.^ We write 

a = aeven + O'CP, (2.3) 

where the non-interesting CP-even part of the cross section (Jeven results from |Alo|^ and 
from the terms in that are independent of 7cp or even in this quantity. The CP- 

violating part of the cross section can be written in the most general Lorentz-invariant 
form as 

acp oc $i[P, • si - (1 ^ 2)] + $2[Pt ■ si - (1 ^ 2)] + ^s[R ■ si - (1 ^ 2)] 
+ *i[(P • si){Pg ■ S2) - (1 ^ 2)] + *2[(P • si)(Pt • S2) - (1 ^ 2)] 
+ *3[(P.-«i)(Pf«2)-(l--2)], (2.4) 

with the auxiliary pseudovector 

Ra = e^^.^Pg^Q'-Pf. (2.5) 

^When calculating linear terms in jcp in we only considered the contribution from diagram h 

in fig. 3, which is important for uih above the if threshold [4]. 
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Figure 3: Feynman diagrams of the underlying gg ti reaction with neutral non-Standard- 
Model Higgs exchanges (dashed). The crossed diagrams are not shown (diagram {g) has 
no crossed partner). 



The coefficients $i, ... $3 and ^^i, . . . ^^3, the structure of which are determined by the 
actual Higgs exchange diagrams (for details, see [3, 17]), depend on s, z and the masses of 
the Higgs boson and the top quark, mn and m^. 

2.2 The Bernreuther— Brandenburg decomposition 

Bernreuther and Brandenburg describe the CP violation in the process of fig. 1 in terms 
of the tt production density matrix, eq. (2.8) in [3] as: 



T^cp = {h^glk + h'g^pi + h'gi'hi) {a' ®1-1®(t') + eijkicgih + Cg2Pi + Cg^hi)a^ ® a\ (2.6) 



uCP^ 



kCP 
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where (in our notation) k = Pg, p = and n is the unit vector in the direction of 
n — PgXPt defined in the gluon-gluon center of mass frame.^ The cr' are Pauh matrices 
with Si = |cr (g) 1 (s2 = |l (8) cr) the spin operator of the top (anti-top) defined in the top 
(anti-top) rest frame. 

The relation of this expansion to our eq. (2.4) is given by: 



s 



s 



\fl j3'A 



-92 



-93 



2mt 



[^*i+/3*2]-s'/3(l--2')*i, 



py/T^^s- (2.7) 



2mt 

The symmetry properties of these functions are given in [3] . 

We confirm the results for these functions (2.7) given in the appendix of ref. [3], up 
to some misprints mentioned in [4]. Also, we have further reduced the results of the box 
diagram to the basic four-point function together with three- and two-point functions. 
This yields a more unique representation, and is also convenient for the use of standard 
loop function libraries, like the "LoopTools" package [13,14]. Our results are collected in 
Appendix A. 

3 Model-independent results 

Before invoking concrete models, we shall first update and review some CP- violating quan- 
tities discussed in [3], both at the parton level and at the hadronic level. 

■^For the momenta of the partons involved, wc use the notation g{ki) + 5(^2) — > t{Pi) + i(P2), while 
in [3] the notation g{pi) + g{p2) t{ki) + i{k2) was used. 



6 



3.1 CP violation at the parton level 

We first consider the non- vanishing^ parton-level CP-odd quantities [3]: 



(P- (Sl - S2))g 
(p- (Si X S2))g 



4 J^^dzAa 
2j'_^dz{zc^f + c^i') 



(3.1) 
(3.2) 



4j^^dzA9 

where is the spin-independent part of the CP-even part of the production density 
matrix given analytically in [3]. We note that these two spin-spin correlations are given 
by absorptive and dispersive parts of the amplitudes (see Appendix A). Plots of these 
expectation values vs. the center of mass energy are shown in fig. 4. 
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Figure 4: Parton-level spin-spin correlations (3.1) and (3.2) in gg tt vs. V^, for 7cp = 1 
and different Higgs masses. Dashed: tjih = 100, 150, . . . , 300 GeV; solid: 350, . . . , 
500 GeV. 



We confirm the results of [3] for (3.1) (fig. 4, left panel), but differ from their results 

for (3.2) (fig. 4, right panel), since there were sign misprints in the structure functions Cg{ 

( f) 

and Cg2 and a wrong factor of Svr in the Higgs width Tz used in [3] and corrected in [4]. 
In our numerical work, we used the "Loop Tools" package [13,14]. 

"'The corresponding quantities involving the gluon or beam direction, vanish due to the symmetry 
properties under z. 
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As stressed in [4], for Higgs masses above the tt threshold, the quantity (3.1) (and to a 
lesser extent (3.2)) has a characteristic "resonance-interference" shape, with peaks in the 
region ~ ttih- As one goes from ^ tjih to \/§ ^ ttih the expectation values (3.1) 
and (3.2) change sign. Thus, there are cancellations which tend to reduce the CP-violating 
effects, when folded with the gluon distribution functions, i.e., when integrated over the tt 
invariant mass, — V§. 

3.2 Observables in pp tiX l+l'X 

We consider events where the t and/or the i quarks decay semileptonically: 

t l+uib, t I'Dih, (3.3) 

and denote the lepton laboratory-frame momenta and energies by l"*^ and E±. Following 
ref. [3], we neglect any CP violation in the top decay. This effect has been found to be 
small [15,16]. 

For reference, we shall consider the observables [3] 

A = E+-E_, (3.4) 

^2 = Pfi'"-Pt-r, (3.5) 

where pt (pf) is the top (anti-top) three-momentum in the laboratory frame. The lepton 
energy asymmetry, Ai, requires events where both the top and the anti-top decay semilep- 
tonically. The observable A2 requires the top or antitop to decay hadronically, while the 
other decays semileptonically. From the two-sample scenario discussed in [4] with sample A 
consisting of events where the t decays semileptonically and the t decays hadronically and 
sample A where the t decays hadronically and the t decays semileptonically, one defines 
the expectation value of the asymmetry A2 as 

(>l2) = (Prn^-(Pfl-)^. (3.6) 



In fig. 5 we show the expectation value (Ai), together with the ratio^ 



(3.7) 



^We shall loosely refer to this as "signal-to-noise ratio" [3] . 



both as functions of the Higgs mass, and for 7cp = 1. The right panel here corresponds to 
fig. 7 of ref. [3], but for nit = 175 GeV.^ Here, and throughout this paper, we have used 
the CTEQ6 parton distribution functions [12]. The fact that S/N = O{10^^) means that of 
the order of a miUion ti events are required (in the dilepton channel). Also, it means that 
the relative lepton energy must be known at (or better than) the per mille level. 
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Figure 5: Left panel: Lepton energy correlations {Ai) in pp —>■ ttX vs. Higgs mass for 
7(7P = 1. Dashed curves labeled "a", "b" and "c" refer to modified observables of eq. (3.8). 
Grey (yellow) band: (A^^^), with 10% uncertainty in Mu- Right panel: Corresponding 
signal-to-noise ratio. 

Next, we turn to the observable A2. We reproduce the analytical result of [3], their 
eq. (4.5), but our numerical result, shown in fig. 6, has the opposite sign, and the signal- 
to-noise ratio is in magnitude smaller than theirs by a factor of (9(1/10) [17]. Thus, the 
observable A2 may be rather hard to access at the LHC 

The reason for the small signal-to-noise ratio for A2 (as compared to Ai) appears to 
be a combination of two effects. The main effect is due to the Lorentz transformations. 
Since it is bilinear in momenta, factors like B2 = 1/(1 — Pgg) and B2 = Pgg/i^ — Pgg) are 
involved, vs. Bi = 1/ ^/l — (3^^ for Ai, with Pgg the gluon-gluon cm. velocity w.r.t. the 
laboratory frame. The "signal-to- noise" ratios for these quantities ( "quasi-observables" ) 

^We reproduce their results for the same choice of mt and the same Higgs widths. Also, we note that 
a factor of (3 is missing in eq. (4.4) of [3]. 
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Figure 6: Left panel: Momentum-momentum correlations {A2) (see eq. (3.5)) in pp ttX 
vs. Higgs mass for 7cp = 1. Right panel: Corresponding signal-to-noise ratio. 



B2 and B2 are reduced with respect to that of Bi by a factor of 3-4. This is due to the 
importance of events where the tt invariant mass is close to threshold, and the relatively 
light it system has a high velocity in the laboratory system. Furthermore, in A2, higher 
powers of cosines of angles are involved, and we note that, for example, (cos^ 6) = 1/3, 
whereas (cos*^^) = 1/5. Together, these two effects give a reduction by (9(1/10). 

For values of the Higgs mass above the tt threshold, there is a fall-off in the absolute 
value of S/N. This is in part due to cancellation between the positive and negative parts 
in (p ■ (si — S2)), as shown in fig. 4. One possible way to enhance the signal, would be 
to consider bins in Mu, centered around the resonant contribution of the (lightest) Higgs 
particle [4]. It would be worthwhile to perform a detailed simulation study of the gains 
of such binning in M^. The value of this technique would depend critically on the mass 
resolution that can be achieved by a given detector. 

We have investigated various other ways to reduce this loss of sensitivity, by studying 
instead of Ai, the modified observables (still for events with two leptons): 

At^^^^^A, (3.8) 
Af ^sigTi{Mu-MH)A^, (3.9) 



10 



= -\Ai\. (3.10) 

The expectation values and signal-to-noise ratios for these observables are also shown 
in fig. 5, labeled "a", "b" and "c". For Af\ and for high Higgs masses, a significant 
improvement is obtained. However, the modified observables A^^^ and Af^ require precise 
knowledge of the tt invariant mass, as well as of the Higgs mass. Of these, presumably the 
tt invariant mass will be the most difficult one to determine. We show in fig. 5 the ranges 
in (Af*) and signal-to-noise ratio that result from errors of 10% in this invariant mass. 

If only the magnitude of Ai is required, then it can be determined with better sensitivity 
from A^-^^ = — l^il, which in fig. 5 coincides with the lower boundary of the shaded region 
and hence gives the largest signal enhancement. 

4 The Two-Higgs-Doublet Model 

In the previous discussion, the amount of CP violation was given by the Yukawa couplings 
a, a (in particular, via the product, 7cp = —aa) and the Higgs mass. In an explicit 
model, the situation is more complex, since there are several Higgs bosons, whose masses 
and couplings will be inter-related by the specific model. As an example of such a more 
complex situation, we shall in this and the next section study a specific Two-Higgs-Doublet 
Model where the CP violation is minimal in structure. 

4.1 Parametrization 

The Two-Higgs-Doublet Model we want to consider, is the one discussed in [6] , where the 
potential is given by 

y = y(0l0if + y(0^02)' + A3(0l0i)(402)+A4(0l02)(40l) (4-1) 



1 

+ 2 



A5((/'l02)' + h.c] - i {m2,(0l0i) + [m?2((/)l(/)2) +h.c.] +m^2(402)} ■ 



The parameters A5 and allowed to be complex, subject to the constraint 

Im = Im A5 V\V2^ (4.2) 

with Vi and V2 the vacuum expectation values (f^ + t>| = t>^, with v = 246 GeV). It is 
this quantity, Immf2 (or ImAs) which leads to CP violation, and one may think of CP 
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violation as being introduced softly, via the mass term 777,^2 in (4.1). Since the potential 
has a Z2 symmetry that is only broken softly by the term, flavour-changing neutral 

currents are suppressed [18,19]. 

The neutral-sector mass squared matrix corresponding to the potential (4.1) can be 
written as (for details, see [6]) 

/ AiC^ + i/s| (A345 - i^)C/3S;3 -ilmAgS^X 



M = v 



(4.3) 



(A345 - v)CjiSii \2S\ + Z/c| -|lm A5 C/3 

\ -ifmAsS^ -iJmAsC^ -ReAs + i// 
with the abbreviations — cos (3, sp — sin/3, tan (3 — V2/V1, A345 = A3 + A4 + ReAs, 
v = Re 777^2/ (2"?^^ /5 cos /3) and /x^ = v'^u. 

For lmA5 ^ 0, the elements Aiu and A423 provide mixing between the 2x2 upper 
left part of the matrix and Aiss- These two sectors would otherwise represent two scalar 
(usually denoted h and H) and one pseudoscalar Higgs boson (denoted ^4). We note that 
these two "mixing elements" are related via tan/3: 

Mi3 = tSLU. P M23- (4.4) 

It is in this sense that the CP violation is minimal in structure. This simple relation is 
violated in more general models, with so-called Xq and A7 terms in the potential [6]. 
We diagonalize (4.3) with the matrix R, defined such that 

RMR^ = diag(Ml^ M^, M|), (4.5) 

and use for the rotation matrix [20] 

cos«6 sinctfeX / cos a sin a 0^ 

1 ~ ^'^^ '^'^'^ ^ 

- sin 0:5 cos ai)j \ Ij 

CaCb S^Cb Sb \ 

^ { -iCaSbSc + SaCc) CaCc-SaSbSc C(, Sc (4.6) 
~l~ ^ce ('^a ~l~ ^c) ^b / 

with Cj = cosftj. Si = sinftj. The angular range, beyond which R is repeated, can be 
chosen as — 7r/2 < a < n/2, — tt < < tt, and —7i/2 < ac < 7r/2. However, the physical 
range is more restricted, since we require Mi < M2 < M3. 
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There are certain symmetries, under which the masses are unchanged, but one or two 
rows of R (i.e., physical Higgs fields) change sign. These symmetries are: 

A : a fixed, ac -(Xc, ai, ^ tt + aj, : Ru ^ -Ru, R2i R2i, Rzi -Rsi, 

B : a ^ TT + a, ab ^ -ab, ac -ac : Ru -Ru, R2i -R2i, Rsi Rsi- 

(4.7) 

In addition, we have the following symmetry for a, ac fixed: 

C : ab ^ TT - ab, ab > {ab ^ -Tr-at, ab < 0) : Ru -Ru, Rjz -Rjz, (4-8) 

the other elements in R being unchanged. Finally, we have, for ac fixed: 

D : (5 — > TT + a : Rji — > -Rji, Rj2 ~^j2, Rjs Rj3- (4.9) 

(The symmetries B and D are of marginal interest, since they only relate the edges of the 
angular range.) 

The CP-conserving case is obtained by taking = or a;, = n, together with etc = 
and a — jTT + a arbitrary. Here, a is the familiar mixing angle of the CP-even sector. 
Thus, as alternatives to ImAs, the angles ab and ac parametrize the mixing that leads to 
CP violation. Of course, replacing one parameter by two, constraints are imposed on the 
mass spectrum. 

Yukawa couplings 

With this notation, and adopting the so-called Model II [21] for the Yukawa couplings, 
where the down-type and up-type quarks are coupled only to 4>i and 4>2, respectively, the 
Hit couplings can be expressed (relative to the SM coupling) as 

Hjtt : -J—r [Rj2 - ij5 cos (3Rjs] = a + iaj^. (4-10) 
sm L/ 

As mentioned in the Introduction, the product of the scalar and pseudoscalar couplings, 

cos (3 / , \ 

7cp = -a a = . 2 o ^J^^js (4-11) 

sm p 

plays an important role in determining the amount of CP violation. We note the following 
symmetries of jcp for any of the three Higgs bosons Hj {j — 1, 2, 3): 
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• under the symmetries "A" and "B" of (4.7), 7cp is invariant, 

• under the symmetries "C" and "D", 7cp changes sign. 

As was seen in sect. 3, unless the Higgs boson is resonant with the tt system, CP 
violation is largest for small Higgs masses. We shall therefore focus on the contributions 
of the lightest Higgs boson. Hi. For the lightest Higgs boson, the coupling (4.10) becomes 



1 r • ~ / \ • n ■ / \i -1 1 sin a sin(2Q;b) 

^— sm a cos(q;;,) - i-f^ cos /^/sm («{,) , with 7cp = 5 — „ . „ ; 

sm p ^ tan p sm p 



where a and «5 are mixing angles of the Higgs mass matrix as defined above. 

When the three neutral Higgs bosons are light, they will all contribute to the CP- 
violating effects. In fact, in the limit of three mass-degenerate Higgs bosons, the CP 
violation will cancel, since [cf. eq. (4.11)] 

E7cP=^E^.-^.-3 = (4.13) 
^ sm p ^ 

due to the orthogonahty of R. 
"Mclximal CP violation" 

The |7cp| of eq. (4.12) is readily maximized w.r.t. the rotation angles a and 



whereas it increases with decreasing values of tan/3. Negative values of etc are related 
to positive values through the symmetry "A" of eq. (4.7) and need not be considered 
separately. Similarly, the case of ccb = 37r/4 is related to that of ctft = 7r/4 by the symmetry 
"C" and just leads to an over-all change of sign for 7cp. We shall somewhat imprecisely 
(since also H2 and contribute) refer to these cases (4.14) as "maximal CP violation". 

More generally, large CP- violating effects are expected when |Im A5I is large, or, when 
the elements M.13 and M.23 are large, see (4.3). These can be expressed in terms of the 
rotation matrix and the mass eigenvalues as 

Mi3 = RiiRi^Ml + R-2iR2zMl + RuRzzMl 
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M23 = RuRiaMf + R22R23MI + R32R33MI (4.15) 

In this description, one measure of "maximal mixing" would be to ignore the lightest 
mass, and then require the coefficients of M| and M| to be equal in magnitude. This 
would require |sin(Q;)| = or |cos(q;)| = 0, i.e., a — 0, ±|7r or tt, and simultaneously 
I sin(2Q;6)| = | sin(2Q;c)| = 1- We shall refer to these cases as "maximal mixing". 

The bounds on the electric dipole moments of the neutron and the electron impose 
restrictions on the allowed magnitude of CP violation in the 2HDM. We note that, within 
a consistent framework, and for maximal CP violation in the Higgs-gauge sector [22] 
(which amounts to maximizing the product of the three couplings (4.26), and thus is 
different from the present concept of "maximal CP violation", where we consider Yukawa 
interactions only), those bounds restrict the mass sphtting between M2 and M3 to be less 
than 0(15% -20%) [23]. 

4.2 Physical content 

Specifying all the parameters of the potential, as well as the structure of the Yukawa 
couplings, the physical content of the model is fixed. We shall here follow a somewhat 
different approach: We start out by specifying the masses of two Higgs bosons (the lightest 
and one other), together with tan/3 and the three angles that determine R. Then, unless 
ab — or ac — 0, the third Higgs mass can be determined, as well as all the couplings. 
This approach gives more "control" of the physical input. 

With tan/3, R and the masses of two neutral Higgs bosons fixed, the third one is given 

by 

RisiRu - R12 tan /3)M^ + i?23(i?2i - R22 tan /3)M| + ^33(^^31 - ^^32 tan /3)M| = 0, (4.16) 

where we have assumed ImAs 7^ and made use of the relation (4.4). Invoking the 
orthogonality of R, one sees that this relation (4.16) only relates differences of masses 
squared: 

i?23(i?2i - i?22tan/3)(M| - M^) + i?33(i?3i - i?32tan/3)(M| - Mf) = 0. (4.17) 

In general, as two Higgs masses approach each other, also the third mass has to approach 
the same value, as we can read off eq. (4.17). However, for particular choices of the 
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parameters, this is not the case. For example, if we put the coefficient of in eq. (4.16) 
equal to zero, one has M2 — M3 for 

Ri3 {Rn- R12 tan /9) = 0. (4.18) 

This equation is satisfied for (i) Ri^ = 0, or (ii) tan/3 = R11/R12 or (in) Ru and R12 each 
zero. Solution (i) implies — 0, whereas solution (ii) implies tan/3 = cot a, or a — |7r — /3. 
Solution(iii) implies \ab\ — 7r/2 where, non- interestingly, the three masses are arbitrarily 
chosen. Following from (ii), for example, we get M2 = M3 for arbitrary Mi, for the choice 
of parameters tan/3 = 1 and d = n/i. 

Allowed regions in the pEtrameter space 

In fig. 7 we show the allowed regions in the a^-ac plane for selected values of tan/3 and a, 
subject to the constraint Mi < M2 < M3. Regions of negative ccc and \ab\ > 7r/2 are not 
shown, they follow by the symmetries "A" and "C", respectively, of eq. (4.7). 

In this analysis, we keep Mi and M2 fixed, and determine M3 from eq. (4.16). Clearly, 
this breaks down when 

i?33(i?3i-i?32tan/3) = 0, (4.19) 

in which case M3 is not determined. This occurs deep inside the unphysical region, as 
explained below. 

The dark (blue) regions in fig. 7 are allowed, the white ones are not allowed. Some 

reasons why a region is physically forbidden are: {{) the equation for M3 has no solution (for 
example, when M| < 0, or when eq. (4.19) holds), {ii) M3 > M2 is violated, {Hi) positivity 
of the potential (4.1) is violated, or {iv) perturbativity is violated. The requirement to 
perturbativity is taken as 

\Xi\ < AttC, (4.20) 

where we somewhat arbitrarily take ^ = 0.8. The lightly shaded (yellow) regions arc also 

allowed, if we use a less stringent condition on perturbativity, ^ = 1.0.'' The cut-off against 

''These numbers are of course only order-of-magnitude indicators. For example, one could argue that 
since Ai and A2 are accompanied by factors 1/2 in the potential (and hence in all couplings), the relevant 
limit is Stt, not An. 
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the disallowed region depends on the additional parameters, which we in this figure take 
as Ml = 100 GeV, M2 = 300 GeV, Mh± = 500 GeV, and /i = (corresponding to Case I, 
studied in sect. 5.1). This dependence enters via the conversion of the masses of the neutral 
Higgs bosons, and the rotation angles, to A's. 



tan/? = 0.5 tan,(?=1.0 tan,i5 = 2.0 




Figure 7: Physically allowed regions in the plane, for tan/3 = 0.5, 1 and 2 and for 

selected values of a. Blue (dark): allowed regions; white: not allowed; yellow: marginally 
allowed (see text). The dashed lines at = ±7r/4 are lines of "maximum CP violation". 
Solid contours in allowed regions: M2/M3 = 0.6 and 0.8 (M2 = M3 along = 0). For no 
choice of mass parameters can the allowed region extend beyond the solid contours in the 
forbidden region. 



The condition (4.19), which is satisfied for 

(1) i?33 = 0, or (2) i?3i-i?32tan/3 = 0, (4.21) 
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defines borders where M| changes sign. We note that these borders are independent of 
the mass parameters of the model. Thus, they are absolute borders which the physically 
allowed region can never cross for any choice of mass parameters. Let us now discuss the 
two cases of eq. (4.21). In the afe-cic plane, case (1) is satisfied when cos a;, = or when 
cosOc = 0. When cosa^ — 0, then eq. (4.16) is trivially satisfied, so this is not interesting. 
The second subcase, cos etc = 0, corresponds to the upper boundary of the plots in fig. 7, 
ccc = 7r/2. Since this line represents a border where M| changes from +oo to — oo, there 
is a region next to it where M3 is so large that perturbativity breaks down (and hence 
forbidden), as indeed seen in the figure. 

Coming now to the more interesting case (2) of (4.21), which, according to eq. (4.6) is 
satisfied for 

tan(/3 + d) tan etc — sin ah, (4.22) 

we obtain the contours shown in fig. 7 deep inside the unphysical region. 

These contours must be in the same "quadrant" in the a^-ac plane as where the phys- 
ically allowed region is located. This can be seen as follows. On one side of the contour, 

M| is large and positive. By continuity, moving from the contour toward the line ab = 0, 
where M3 = M2, the value of M3 will inevitably pass through an allowed region where 
Ms Z 0{M2). 

This observation allows us to determine in which quadrant the physically allowed regions 
are located. Near the origin, (4.22) can be approximated as 

tan(/3 + a)~— . (4.23) 

Thus, abOLc > when 

< /9 + a < Itt, First (or third) quadrant (4.24) 

whereas abO-c < when 

— ^7r<P + d<0 or < P + d < Second (or fourth) quadrant (4-25) 

in agreement with fig. 7. Finally, we note from (4.23) that the physical region may occupy 
a large fraction of the actual quadrant when |/? + cil ~ 7r/2, and only a small fraction of 
the quadrant when |/? + q;| ~ or tt. 
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Higgs— vector-vector couplings 

The non-discovery of a Higgs boson at LEP poses constraints on tan f3 and on the Higgs 
mass. While SM Higgs masses below 115 GeV are excluded [24], as well as tan/? ^ 2-3 
in the MSSM [25], these bounds can be eluded in the 2HDM, if the hghtest Higgs boson 
couples with sufficiently reduced strength to the vector bosons. In fact, this coupling is 
(relative to the corresponding SM coupling) given by 



H,ZZ 



cos P Ril + sin P Ri2. 



(4.26) 



Higgs masses down to about 50 GeV are allowed, provided this couphng squared is less 
than about 0.5 [26]. 
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Figure 8: Squared coupling of the lightest Higgs boson to vector bosons (relative to that 
of the Standard Model), Qyy, vs. mixing angle a^, for fixed a = 0, 7r/4 and 7r/2. These 
couplings are symmetric under 7i — aii, and independent of etc- For a = and tan/? = 
0.5, only small values of at are allowed for the chosen mass parameters: Mi = 100 GeV, 
Ms = 300 GeV, Mh± = 500 GeV, = tt/S and ^ = 0.8 (cf. fig. 7). 



We show in fig. 8 the square of the coupling of the lightest Higgs boson to vector bosons, 
relative to that of the Standard Model, Qyy, for low values of tan/3. As discussed above, 
for masses of the lightest Higgs below ~ 100 GeV, this quantity should be well below unity, 
in order not to be in conflict with the non-observation of a Higgs boson at LEP2 [24]. 

It should be noted that this coupling can be expressed in terms of angles only, without 
any dependence on masses. However, for certain combinations of mass parameters and etc 
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values, some ab ranges might be unphysical, as illustrated for o; = and tan/3 = 0.5 in 
this figure, as well as in fig. 7. 

5 Case Studies 

In order to get a better idea how much CP violation the 2HDM can give, we shall here 
discuss the case of "maximal CP violation" in the sense of eq. (4.14), together with small 
values of tan (3 and light Higgs masses. 

In the following two subsections, we shall analyze two cases, first (sec. 5.1) the case 
of two Higgs bosons being light, and then (sec. 5.2) the case of one fight and two heavy 
ones. These are quafitatively different, since in the first case there will be considerable 
cancellations among the CP- violating contributions, because of (4.13). 

With all mixing angles fixed, and two Higgs masses specified, the third Higgs mass is 
determined. Thus, the "soft" mass term /i^ = Re 771^2/ sin 2/3 only enters in converting 
the masses (and angles) to A's. Wc check that the required A's satisfy positivity of the 
potential and perturbativity, eq. (4.20). 

5.1 Case 1. Two light and one heavy Higgs bosons 

As a first case, we consider Mi < M2 — 300 GeV (below the tt threshold) with Mh± — 
500 GeV and tan /3 = 0.5 and 1.0. (This is well within the fimits on Mh± and tan (3 derived 
from studies of meson decays and mixings [27].) Then, for ac = vr/4, 7r/3 and 57r/12, we 
show in fig. 9 the heaviest Higgs boson mass, M3 vs. Mi, for tan/3 = 0.5, 1.0 and 1.5. We 
note that, for some of the parameters considered, there is no solution (we take ^ = 1.0), as 
also follows from fig. 7. 

In fig. 10 we show, for this case of "maximal CP violation" and low values of tan/3 
(tan/3 = 0.5 and 1.0), the corresponding signal-to- noise ratio of the observable Ai (3.4), 
for Ml ranging from 80 to 300 GeV, keeping M2 = 300 GeV fixed and letting M3 vary 
accordingly. We show separately the contributions of only the lightest Higgs boson, and the 
contributions of all three. Considering only the contribution of the lightest Higgs boson, the 
CP- violating effects are seen to be significantly larger than in the "general case" studied 
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Figure 9: Heaviest Higgs-boson mass, M3 vs. Mi, for M2 = 300 GeV, and Mh± = 500 GeV, 
for three values of ac = 7r/4, n/S and 57r/12, and three values of tan/3 = 0.5, 1.0 and 1.5. 
Furthermore, we consider "maximal CP violation" with a = 7r/2 and at = — 7r/4. 

in sect. 3, essentially because of the enhancement of 7cp due to the small values of tan (3. 
We note that, for the parameters considered here {a = 7t/2 and ab = — 7r/4), and for the 
lightest Higgs boson, ■jcp is negative, contrary to the case studied in sect. 3. 

However, in the present case, the inclusion of all three contributions leads to large 
cancellations. This also implies that, even if no CP violation should be observed, it may 
be difficult to conclude that the Higgs sector conserves CP, since such cancellations are 
possible. 

As expected, the CP- violating effects are largest for low values of Mi. Since the two 
lightest Higgs bosons are below the tt resonance, there is a smooth dependence on Mi. 
As Ml approaches M2 (where all three masses are degenerate), the CP- violating effects 
cancel. On the whole, the resulting CP violation is reduced by a factor of ~ 1/4 compared 
to the "model-independent" case considered in sect. 3. 

Whereas the top Yukawa couplings of the lightest Higgs boson do not depend on ac, 
some dependence on this quantity enters, via the couplings of the two heavier Higgs bosons 
to the t quark. This is also illustrated in fig. 10. 

Since the two heavier Higgs bosons can contribute significantly to {Ai), the dependence 
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M, [GeV] 

Figure 10: Signal-to-noise ratio, for lepton energy correlations (Ai) in pp — ^ tiX vs. lightest 
Higgs mass for the 2HDM, for a case of three moderately light Higgs bosons (see fig. 9) 
and tan/3 = 0.5 and 1.0. "Maximal CP violation" is considered (cf. eq. (4.14)), and 
Oc = vr/4, tt/3 and 57r/12. Also shown, is the contribution of the lightest Higgs boson only 
(independent of etc). 

on ttfe need not be as simple as discussed in sect. 4.1. In particular, the CP violation need 
not be maximal for af, = ±7r/4 (or ±37r/4). Fig. 11 (left panel) illustrates the dependence 
on ab for fixed a = 7r/2, ac = vr/3, and two values of tan/3 (0.5 and 1.0). One confirms 
the general impression from the form of (4.6) and the discussion following (4.15), that \ab\ 
should be sizable, but not necessarily tt/A, in order to maximize the CP violation. 

It is also possible that one Higgs boson does not violate CP in its couplings to the t 
quark, and yet the other two do. For example, for a = 0, the lightest Higgs boson has 
a pure pseudoscalar coupling to the top quark, and hence the CP violation is exclusively 
due to the two heavier Higgs bosons. In this case, the resulting total contribution to the 
lepton energy correlations (Ai) may even exceed that obtained when the contribution of 
the lightest Higgs boson is maximal, see fig. 11 (right panel). 

5.2 Case 2. One light and two heavy neutral Higgs bosons 

Next, we let two of the neutral Higgs boson masses be more heavy, taking M3 > M2 = 
500 GeV, and M{j± = 700 GeV. The resulting mass values, M3, are shown in fig. 12, for 
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Figure 11: Sensitivity of Ai in pp — > tiX vs. lightest Higgs mass. Left panel: Different 
values of ab are studied, for fixed a = 7r/2, ac = n/S. Right panel: Different values of a are 
studied, a = 0, 7r/4 and tt/2 ("maximal mixing") with etc = 7r/3 (see fig. 7). The mixing 
angle \ab\ — 7r/4 is kept fixed. Two values of tan/3 are considered: 0.5 and 1.0. 



two values of a, three values of etc = 7r/4, tt/S and 57r/12, and small values of tan/3. For 
the cases studied in sect. 5.1, solutions could be found with n — 0. This is here not the 
case. For a — 7r/2 ("maximal CP violation"), no value of // allows a hght Hi (fig. 12, right 
panel), so we have included also a — ti/3 (left panel), where light Hi can be realized for 
large /i. 

We show in fig. 13 the corresponding results for the lepton energy asymmetry, (^i), 
for tan/3 — 0.5 and 1. In this case, as opposed to the cases considered in sect. 5.1, since 
H2 and H3 are heavier, it is a good approximation to consider only Hi exchange up to 
Ml ~ 450 GeV. As a result, the CP- violating effects are significantly larger than in the 
low-mass case studied in sect. 5.1. The amount of CP- violation for one light Higgs boson. 
Ml = 0(100 GeV), is comparable with, and may exceed that obtained for a mass around 
the ti resonance. 

We note that there is only a rather weak dependence on ac- This is clearly due to the 
fact that the contribution of Hi (which is independent of etc) dominates. (However, as 
etc — > 0, the physical range of Mi shrinks, cf. fig. 7.) 
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Figure 12: Heaviest Higgs-boson mass, M3 vs. Mi, for M2 = 500 GeV and Mh± = 700 GeV. 
For some parameters, there is a lower limit for Mi, due to perturbativity, eq. (4.20), where 
we take ^ = 1.0. 



Fig. 14 is devoted to a study of the dependence of ab- As expected (since Hi dominates), 
the CP violation is maximized for at — —it /A. As ab ~^ 0, the effect vanishes. There are in 
fact two reasons for this. First, the contribution of Hi to 7cp vanishes linearly. Secondly, 
M2 and M3 become degenerate in this limit. According to the discussion in sect. 4.2, all 
three masses must then become degenerate, and there is no CP violation. 

5.3 Case 3. The transition region 

As is seen from fig. 4, and also from a qualitative comparison of Case 1 and Case 2, 
there is in general less CP- violation if two Higgs masses are below the ti threshold than 
if two masses are above. We here explore this transition region in a little more detail, 
by comparing in fig. 15 a few values of the intermediate Higgs boson mass M2, below and 
above the tt threshold, keeping the mixing angles fixed at a = 7r/2, ab = — vr/4, ac = vr/3 
and tan/3 = 0.5. One sees that, for values of Mi well below M2, there is little dependence 
on M2, unless M2 is above the tt threshold, in which case the CP-violation (for fixed Mi) 
increases rapidly with M2. 
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M, [GeV] M, [GeV] 

Figure 13: Signal-to-noise ratio, for lepton energy correlations (Ai) in pp tiX vs. Higgs 
mass for the 2HDM. Similar to fig. 10 for a = n/3 and M3 > M2 = 500 GeV, with 
Mh± = 700 GeV. Note that one set of curves is "upside-down" with respect to the other, 
since the aft-values have opposite signs (we keep etc positive in both cases). 



6 Concluding remarks 

We have reviewed and confirmed the basic results of [3] on CP violation induced by non- 
standard Higgs exchange in the production of ti pairs at the LHC, up to a few misprints. 
The results for the box diagram have been further reduced to basic loop integrals. Some 
modifications of their observable Ai = — (where E± are the lepton energies) have 
been investigated. These modifications can in some kinematical situations improve the 
signal-to-noise ratio, provided the ti invariant mass can be determined with a reasonable 
precision. 

For the observable ^2 = pf • 1^ — ■ I~ , we find a much lower signal-to- noise ratio than 
was given in [3]. This observable thus appears less suited for a search for CP violation at 
the LHC (for details, see sect. 3.2). 

In a minimal, CP- violating Two-Higgs-Doublet Model, where couplings and masses are 
constrained and related by the Higgs potential and the Yukawa couplings ("Model II"), the 
CP-violating effects can at the LHC be at the per mille level, as given by Bernreuther and 
Brandenburg [3] (see also [1]), provided that (i) tan/3 is small (^ 1), and that (ii) there is 
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Figure 14: Signal-to- noise ratio, for lepton energy correlations (Ai) in pp tiX vs. Higgs 
mass for the 2HDM. Similar to fig. 13 for three values of ab- 
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Figure 15: Lepton energy correlations (Ai) in pp —>■ ttX vs. Higgs mass for the 2HDM. 
The intermediate Higgs mass, M2, is in the tt threshold region. 

only one light Higgs boson. If there are two neutral Higgs bosons below the ti threshold, 
the effects are severely reduced by cancellations among different contributions. 

Our study does not include QCD corrections. Within pure QCD (no Higgs exchange), 
next-to-leading-order corrections are known [28]. They are not Isnown for these CP- 
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violating effects. Since we study ratios, one may hope that the dominant parts of the 
QCD corrections cancel, but this remains to be seen. 

Finally, we recall that there are other observables that we have not studied (see, e.g. [4, 
29]). At the parton level, the basic asymmetries are larger, of the order of 10~^, as compared 
with the order 10~^ effects discussed here. Thus, for some other choice of observable, it is 
quite possible that a larger fraction of the 2HDM parameter space could be explored. 
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Appendix A 

In this appendix we collect complete results for the functions fe^/', ■> Cgi, Cg2 of cq. (2.6), 
correcting some misprints in [3] and also reducing completely the contribution from the box 
diagram to four-, three-, and two-point scalar integrals. We define the following two-point 
functions: 

(A.l) 

Bo{s,mlmf) = / A r._. , i^x, IlaTFT: TT^o 1727' (^-2) 
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[(? - iQr ■ 


-mj\[{q- 
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-mj\[{q- 


- iPt)' - ml] 
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three-point functions: 

If 1 



ln\2 _ ' 



m 

If 1 

(A.5) 
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If 1 

C,{t,m%,mlm',) = J d\ _ , ^^^^ _ _ , _ _ , ^^^^ _ , 

(A.6) 

and the basic four-point function: 

I [ n 1 

^oW - y '^"^[(g + \QY - ml][{q - \P,Y - ml][{q - \QY - m?][(g - \P,)^ - m\\ ' 

(A.7) 

Here, s — Q'^, t — —Q^{1 — Pz). Furthermore, for u — —Q^{1 + Pz), the scalar integrals 

BQ^u^m^j^ml), Co{u,rn'jj, 171^,171^) and Dq{u) are obtained from (A. 3), (A.6) and (A.7), 
respectively, with Pg — * —Pg (or Pf —Pt)- Note that, in this notation of [3], the function 
Co{s,mf,mf,mf) can not be obtained from Co(s, m^, m^, m|) by a simple substitution. 
We also define the abbreviations 

Dpii) = mjjDo{i) - (1 - fiz)Co{i,ml,mlm''t), 

Df{u) = m%Do{u) - (1 + pz)Co{u, m^, m^^ m^^), (A.8) 
where /5 = (1 — 4m|/s)^/2 and z = {Pg ■ Pt), together with the UV- finite combination 

B = 9Pz[BQ{i,mH,mf) - Bo{u,mjj,mf)] + 18z^[Bo{s,m'^,m^) - Bo{m^,rnjj,rn^)] 
+ 7[Bo(t, mjj, ml) + Bo(m, m^, m^) - 25o(m2, m|^, m^)], (A.9) 



and (see 



2 2 2N , BQ{ml,m]j,ml) - Bo{i,m%,m 



CM = Co(t, m^, m^, m^) + " pv-.-^.-t. ^ (A.IO) 



ml — t 



-[m],Co{s, m\, ml m\) + i?o(5, m^, m\) - B,{ml m\, m: 



Furthermore, Cs{u) is obtained from eq. (A. 10) by the replacement i ^ u. The factor 

(A.12) 

where gg is the QCD coupling constant, determines the over-all scale. 

The functions b'^-f and Cgi are odd under (z —z, implying i ^ u) whereas 
and Cg2 are even. In the following, the four- and three-point functions are finite in four 
dimensions, as well as all combinations of two-point functions. 
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For the box diagram (diagram (c) in fig. 3), we find: 



-rrit E 



192(1 -/32^2)(i_ ^2)5^ 



X 



(7 + ^(3z) (2z[ s (1 - (3zf + Am]j\ IrnDpii) 

+ {s(3z[s/3{l - (3z) (1 - 2(3z + z^) + 4m^(/? - z)] 
-s^/3(l-/3^)(l-/3^)}ImDo(^)) -{z^-z] i ^ u) 
+ Ap'^z [s{y - 18/5^^2 + 7z^) - 18m^z2(3 - z^) - 9s(l - f)] Im Co(s, m^, m], ml) 
+ Az[s{7-ll(5'^z^) + 2Sm]j\lmCQ{s,ml,ml,ml)-A(5'^z{l-z^)lmB , (A.13) 

, (c) _ 

X [(7 + 9pz) (2[ s (1 - pz)^ + Amjj] [(1 - z'^)mt + E^z^] Im L>F(i) 
+ {s(3[s(3{l -(3z){l- 2(3z + z"") + 4m^(/3 - ^)] [(1 - z'')mt + 

- fEipzil - Pz){l - ImDoit)) + {z^ -z- i u) 

+ 4/32{[s(7 - 18/3 V + 7^2) - 18m|^2(3 - z^)][{l - z^)mt + E^z^] 

- 18z^mt[m|(l - + 2mt£;l]}ImCo(s,m^,mJ^m^) (A.14) 
+ 4[s(7 - 11/3^^2) + 28m^][(l - z^)mt + Eiz'^]lmCQ{s,ml,ml,ml) 

- 4/3^(1 - [(2 - z^)mt + ^i^2]Im^l , 



(c) _ -^1 

Cg^ - K^CP ^g2(l _ /32^2)(1 _ ^2)^-^ 



X 



(7 + %(3z){2z[ s {0^z^ + 2/3^ - 2/3z - 1) + 4m^] Re Dpit) 
+ {s/3z[s/3(l - pz){l - 2(3z + z^) + 4m|(/3 - z)] 

+ 5^/3(1 - /3;2)(1 - /32)} ReDo(t)} - {z ^ -z; i u) (A.15) 
+ Ap'^z [s{7 - 18/3^^2 + 7z'^) - 18m%z'^{3 - z^) + 9s(l - (5'^)] ReCo(s, m|, ml, ml) 
- Az[s{l + 11/3^^2) - 28m^ - 14s/32]ReCo(s, ml, ml, ml) - 4/3^z(l - ^^)ReB] , 

(c) _ . 

192(l-/32z2)(i_^2)^^ 

(7 + 9/3;z)(2{[s(/32^2- 2/3^-1) + Am]j\[{l-z^)mt + + 2sEip'^z^}ReDF{t) 
+ {s/3[s/3(l - (5z){l - 2(5z + + 4m^(/3 - ^)][(1 - z')mt + ^i^'] 



X 
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+ s'^EiPz{l - Pz){l - P^)} ReDo{i)) + {z ^ -z; i u) 
+ 4/5^{[s(7 - 18/3^^2 + 7^2) - 18774^2(3 - z'^)][{l - z^)mt + E^z^] 

- \^z^mt{m\{\ - z^) - 2mtEi]}ReCo{s,7n'jj,m'^t,m'^t) (A.16) 

- 4{[s(7 + nf3h'^) - 28m|][(l - z'^)mt + E^z'^\ - UsE^0^z'^}Re Co(s, m^) 

- 4/3^(1 - z^)[{2 - z^)mt + Eiz^]ReB. 

Prom the vertex diagrams {d and e in fig. 3): 
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(A.17) 
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Prom the self-energy diagram (after fixing the signs): 
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(7 + 9pz) 



El + mt Ef{l - Pz) 



X ^-Bo(m^, mjf, m]) - Bo{i, mjj, m^)^ + {z ^ -z; u) 
Prom diagram (g): 

-3m2/33^2jjj^G'(s) 



(A.18) 



4? = ^CP 



8(1 - /32^2) 

3m|^VReG(£) 

8(1 - /32^2) ' 



(A.19) 



Pinally, from diagram {h): 



hf^^^ilcp 



-P 



{s -m^y + Tlml 4(1-/3 



(A.20) 



X 



2m( [rijmijReCo(s, m^,m^,m^) — {s — m^)ImCo(s, mj] + TuTaH 
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{s - m\Y + V\m\ Ll6(l - /J^^^) 
X [(s — m|f)Re Co(s, m^, m^, m^) + Vjiuiuhn. Co{s, m^, m^, m^)] — 4(s — m^)^ 

3 m^v^ , 
32 Sn^ 



(A.21) 



+ 



where Th is the width of the Higgs boson calculated by summing the partial widths for 
H W~^W~, ZZ, tt in the version of 2HDM we considered. 
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